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Nomenclature 

Dynamic  Mode  Decomposition 
Gas-Centered  Swirl-Coaxial  Injector 
Momentum  Flux  Ratio 
Pressure  Node 

Proper  Orthogonal  Decomposition 


I.  Introduction 

Gas-centered  swirl-coaxial  injectors  (GCSC)  have  seen  an  increase  in  interest  over  the  past  decade  for  rocket  engine 
applications1"4.  These  injectors  combine  features  of  the  coaxial  air-blast  and  pressure  swirl  atomizers.  Similar  to 
coaxial  air  blast  atomizers,  a  shear  layer  is  formed  between  the  gaseous  and  liquid  propellants.  For  these  injectors,  the 
liquid  propellant  is  injected  tangentially  through  the  outer  annulus,  imparting  a  swirl  onto  the  liquid  similar  to  the 
pressure  swirl  atomizers.  The  resulting  shear  and  tangential  forces  are  responsible  for  atomizing  the  liquid  propellant. 
This  injector  is  ideally  suited  for  oxygen-rich  staged  combustion  cycles  and  has  been  used  successfully  in  the  Russian 
RD-1805  rocket  engines.  Although,  great  operational  success  has  been  achieved  with  these  injectors,  there  is  very  little 
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understanding  regarding  the  full  modes  of  operation  these  injectors  experience  at  rocket  relevant  conditions  and  their 
susceptibility  to  combustion  instabilities. 

Most  of  the  related  studies  performed  using  GCSC  injectors  have  been  focused  on  the  atomization  process  of  the 
liquid  fuel  sheet,  hydrodynamic  instabilities,  and  spray  characteristics.  Researchers  at  the  Air  Force  Research 
Laboratory  (AFRL)6"8  have  developed  physics-based  scaling  laws  for  GCSC  injector  design  through  extensive  cold 
flow  experimental  studies  with  a  parametric  variation  of  physical  dimensions  such  as  cup  recess  length  (Lc  shown  in 
Figure  1).  It  has  been  shown  that  the  liquid  sheet  length  prior  to  break-up  scales  strongly  with  the  momentum  flux 
ratio  (J)  of  the  two  propellants.  Lightfoot  et  al.6' 8  showed  that  the  liquid  sheet/film  break  up  length  reaches  an 
asymptomatic  value  for  a  critical  J  value,  hence  any  operating  J  value  greater  than  the  critical  will  result  in  the  same 
film  break  up  length.  Experimental  studies  by  Kim  et  al.9  on  the  spray  patterns  for  a  GCSC  injector  performed  at  high 
pressures  also  identified  J  as  the  controlling  parameter  for  the  break  up  length.  Schumaker  et.al.2  and  Jeon  et  al.10 
studied  the  effect  of  the  recess  cup  length  Lc  for  GCSC  injectors.  They  independently  showed  experimentally  that  Lc 
has  very  little  effect  on  the  film  break-up  length  once  a  critical  operating  J  value  is  reached.  Eberhart  et  al.n 
investigated  the  the  spray  cone  angle  as  a  function  of  chamber  pressure,  observing  that  the  cone  angle  from  the  swirling 
liquid  sheet  decreased  with  increasing  chamber  pressure. 

As  the  bulk  of  the  prior  research  has  been  focused  on  characterizing  the  spray  pattern  and  design  methodologies  for 
GCSC  injectors  under  cold  flow  conditions,  there  is  very  little  understanding  of  how  these  injectors  operate  under 
reacting  conditions.  A  numerical  study  by  Muss  et  al.12  demonstrated  phenomenological  differences  between  cold 
flow  and  hot-fire  spray  patterns,  particularly  droplet  size.  Pomeroy  et  al.13  and  Morgan  et  al.14  studied  the  response  of 
GCSC  injectors  to  transverse  instabilities,  but  no  effort  to  characterize  the  flame  in  steady-state  conditions  was 
attempted  as  the  combustion  chamber  produced  a  self-excited  combustion  instability. 

This  paper  encompasses  an  experimental  effort  to  describe  the  steady-state  combustion  behavior  and  the  combustion 
response  of  a  single  GCSC  injector  under  acoustic  forcing  conditions  using  high-speed  imaging  techniques.  Varying 
momentum  flux  ratios  were  investigated  for  a  single  geometry  GCSC  injector  to  characterize  steady-state  operation, 
whereas  one  test  condition  was  subjected  to  high-frequency,  low  amplitude  transverse  acoustic  oscillations  in  a 
pressure  node  (PN)  configuration  for  the  forcing  study.  All  experimental  efforts  were  done  at  a  high  chamber  pressure 
of  3.2  MPa,  where  simultaneous  high-speed  shadowgraph  and  OH*/CH*  chemiluminescence  images  were  taken  to 
capture  the  flame  behavior. 


II.  Experimental  Setup 

The  GCSC  injector  was  tested  in  the  Combustion  Stability  Lab  (EC-4)  at  Edwards  AFRL  CA,  using  RP-2  and 
gaseous  oxygen.  A  detailed  review  of  the  complete  system  can  be  found  in  the  work  of  Wegener  et  al.15.  The  system 
has  since  been  upgraded  to  accommodate  liquid  hydrocarbon  propellants.  This  test  facility  also  has  fuel  systems 
capable  of  flowing  gaseous  fuels  such  as  CEL  and  EL  and/or  cryogenic  oxygen  or  gaseous  oxygen.  The  facility  is 
equipped  with  a  Pacific  Instruments  data  acquisition  system  that  is  used  for  diagnostic  and  system  monitoring 
measurements.  The  GCSC  propellants  were  ignited  using  a  hydrogen/air  torch  at  the  operating  chamber  pressure  with 
a  capsule  of  nano-particulate  additives  photo-ignited  by  a  Xenon  flashbulb  acting  as  a  pre-ignition  source  for  the  torch. 
A  detailed  description  for  the  ignition  system  can  be  found  in  the  work  of  Badakshan  et  al.16 

The  experiment  was  performed  within  a  high  pressure,  optically  accessible  combustor.  A  CAD  rending  of  the 
windowed  combustor  is  shown  in  Figure  2.  This  windowed  combustor  features  a  chamber  within  a  chamber  design, 
where  the  large  outer  windows  are  used  to  contain  pressure  and  the  small  inner  windows  are  used  to  contain  the  flame. 
The  windowed  combustor  features  two  piezo-sirens  positioned  on  either  side  of  the  injector  that  are  capable  of 
generating  an  acoustic  standing  wave  to  force  the  flame/injector  flow.  The  combustor’s  two  large  windows  are  used 
for  the  various  line-of-sight  imaging  techniques  (e.g.,  shadowgraphy).  The  combustor  also  features  off-axis  windows 
for  future  laser  diagnostics  measurements  like  particle  image  velocimetry  (PIV)  and  planar  laser  induced  fluorescence 
(PLIF).  The  windowed  chamber  is  pressurized  using  a  continuous  flowrate  of  nitrogen  gas  and  is  not  affected  by  the 
combustion  process,  as  the  volume  occupied  by  the  flame  is  small  relative  to  the  volume  of  the  pressurizing  nitrogen 
gas.  This  helps  sever  the  thermo -acoustic  feedback  loop,  effectively  isolating  the  flame  from  any  potential  chamber 
coupled  combustion  instabilities. 

Three  high-speed  cameras,  (Phantom  v710,  vl210,  and  v2512)  were  used  for  the  simultaneous  shadowgraph,  OH*, 
and  CH*  chemiluminescence..  All  three  cameras  were  synchronized  using  an  IRIG  time  signal,  along  with  the  system 
measurements  from  the  Pacific  Instruments  data  acquisition  system.  A  schematic  of  the  optical  setup  for  the  three 
high-speed  imaging  techniques  is  shown  in  Figure  3.  For  the  shadowgraph  system,  a  broadband  collimated  light  travels 
through  two  dichroic  mirrors.  This  collimated  light  first  encounters  a  dichroic  mirror  coated  to  reflect  410  nm  to  450 
nm  light,  and  transmit  in  the  higher  visible  wavelengths,  towards  the  Phantom  v2512  for  CH*  chemiluminescence. 
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The  CH*  chemiluminescence  431  nm  light  is  captured  by  a  Lambert  Instruments  HiCATT  intensifier  with  a  Semrock 
filter  (FFO 1-425/26).  The  reflective  surface  of  this  dichroic  mirror  rejected  the  blue  light  portion  from  the  broadband 
collimated  light  away  from  the  intensifier  and  hence  there  was  no  need  to  pre -color  filter  this  light  to  avoid  saturating 
the  intensifier.  The  second  dichroic  mirror  is  coated  to  reflect  the  UV  light  for  the  OH*  chemiluminescence  and 
transmits  visible  light  used  for  the  shadowgraphy  system.  The  OH*  chemiluminescence  was  also  imaged  using  a 
HiCATT  intensifier  with  a  Semrock  filter  (FFO  1-320/40).  The  shadowgraph  camera  was  set  to  a  gate  of  7  ps,  with  the 
intensifiers  gated  at  110  ns  for  OH*  and  500  ns  for  CH*.  The  cameras  captured  images  at  1000  Hz  for  a  one  second 
of  steady-state  combustion.  For  the  acoustically-forced  conditions,  the  image  capture  rate  of  the  shadowgraph  and 
CH*  systems  was  increased  to  25  kHz  for  0.2  seconds,  while  the  OH*  system  was  set  to  a  lower  frame  rate  to  capture 
one  second  of  the  acoustically-forced  flame. 

As  can  be  seen  in  Table  1,  four  test  conditions  were  investigated  with  different  J  values  of  oxidizer-to-fuel 
momentum  flux.  For  each  test  condition,  the  RP-2  flowrate  was  kept  constant  while  the  gaseous  oxygen  flow  was 
varied  to  achieve  different  J  values.  The  flowrates  of  RP-2  and  gaseous  oxygen  were  measured  using  a  Siemens 
Coriolis  mass  flow  meter  (SITRANS  FC)  and  Porter  Thermal  mass  flow  meter  (Type  123),  respectively.  The 
windowed  combustor  was  pressurized  to  3.2  MPa  using  gaseous  nitrogen  prior  to  ignition.  From  these  four  conditions, 
Test  4  was  selected  to  be  acoustically  forced.  During  these  experiments,  the  piezo-sirens  were  operated  180°  out-of- 
phase  to  form  a  standing  wave  pressure  node  (PN)  condition  in  the  vicinity  of  the  reacting  flow  field,  where  a  PN  and 
corresponding  velocity  anti-node  (VAN)  was  located  at  the  injector  location.  Three  differential  pressure  transducers 
(Kulite  XCL-100)  were  mounted  to  one  of  the  inner  windows  to  measure  the  pressure  amplitude,  phase,  and  shape  of 
the  standing  wave  during  testing.  These  inner  windows  were  also  coated  with  an  oleophilic  coating  that  spread  RP 
droplets  optical  thin,  thus  not  obscuring  the  view  to  the  flame. 

Table  1.  Operating  Conditions  and  Average  Cone  Angles 


Test 

rhRP-2(.q/s ) 

rhnx(g/s ) 

j 

Pc  (MPa) 

Left  Side 
STD  0 

Right  Side 
STD  0 

001 

1.36 

5 

0.92 

5.5 

3.26 

5.27 

10.45 

002 

1.40 

5.67 

0.84 

6.45 

3.3 

5.79 

9.94 

003 

1.40 

4.56 

1.04 

4.4 

3.21 

6.84 

9.90 

004 

1.37 

4.62 

1 

4.88 

3.17 

5.58 

7.88 

III.  Results 


A.  Average  Behavior 

For  each  test  condition,  one  second  of  steady-state  flame  operation  was  averaged  using  each  optical  method  and 
can  been  seen  in  Figure  4.  The  vertical  green  line  shown  in  this  figure  represents  10  mm  of  the  flow  field  beginning 
at  the  nozzle  exit  plane.  For  the  CH*  averages,  the  view  is  partially  blocked  downstream  by  the  differential  pressure 
transducer  used  to  characterize  the  standing  wave  for  acoustic  forcing  studies.  From  both  the  CH*  and  OH*  average 
images,  it  is  observed  that  the  flame  is  separated  from  the  nozzle  for  each  test  condition.  Great  care  must  be  taken  in 
interpreting  this  behavior,  since  the  average  shadowgraph  images  show  the  flame  is  attached  at  the  nozzle  exit.  By 
overlapping  the  CH*  and  shadowgraph  average  images,  as  shown  in  Figure  5,  an  apparent  spatial  offset  exists  between 
where  the  CH*  is  first  observed  and  the  flame  at  the  injector  exit  present  in  the  shadowgraph  average.  Related 
experiments  done  at  AFRL  show  strong  indications  of  burning  in  the  cup  recess.  Those  related  experiments  have 
shown  decomposed  fuel  deposits  on  the  injector  face,  significant  heat  flux  into  the  injector  and  evolved  injector  face 
temperatures,  which  are  all  experimental  findings  indicating  burning  in  this  portion  of  the  injector.  Purdue  University 
has  also  measured  OH*  and  CH*  directly  in  the  injector.  It  is  currently  believed  that  the  flame  is  anchored  within  the 
cup  recess  of  the  injector  for  these  studies,  and  it  is  observed  that  combustion  takes  place  in  the  near  injector  region 
and  at  times  the  chemiluminescence  is  obscured  from  view.  It  will  be  shown  later  from  the  instantaneous  shadowgraph 
images  that  highly  luminous  pockets  of  reacting  RP  are  observed  to  be  attached  to  the  nozzle  exit,  then  are  later 
obstructed  from  view  by  either  a  soot  cloud  or  by  the  liquid  RP  sheet  emanating  from  the  nozzle. 

From  Figure  4,  Tests  2  and  3  both  show  an  asymmetric  flame  structure  present  in  the  average  OH*  and  CH* 
images.  This  asymmetry  is  not  fully  understood  since  a  consistent  asymmetry  is  not  present  for  the  other  two  test 
conditions.  In  the  X-ray  radiography  cold  flow  studies  performed  by  Schumaker  et  al .13,  an  asymmetry  of  the  mass 
distribution  in  the  swirling  spray  was  observed  at  different  operating  conditions  for  a  single  injector.  It  is  possible  that 
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a  similar  mechanism  is  occurring  during  Tests  2  and  3,  where  greater  concentrations  of  fuel  are  present  at  different 
spray  locations,  resulting  in  higher  CH*  and  OH*  concentrations  at  those  locations. 

There  is  also  an  asymmetry  present  in  the  cone  angle  formed  by  the  spray.  In  Figure  6,  the  standard  deviation  of 
the  shadowgraph  images  were  used  to  measure  the  cone  formed  by  the  reacting  flow  field.  The  standard  deviation 
images  were  chosen  over  the  averages  for  the  cone  angle  O  measurement  as  there  were  more  pronounced  edge 
boundaries  of  the  reacting  flow  field  present  in  the  standard  deviation  images.  These  cone  angle  measurements  were 
performed  by  visual  inspection  and  are  presented  in  Table  1.  As  the  cone  angles  significantly  differ  between  the  left 
and  right  spray  edge  boundaries  for  all  test  conditions,  an  unequal  cone  angle  contributing  to  this  asymmetry  is 
observed.  Further  investigation  will  be  necessary  to  address  the  physical  mechanisms  leading  to  this  asymmetry. 


B.  Instantaneous  Behavior 

A  time  series  of  the  instantaneous  shadowgraph,  CH*,  and  OH*  images  for  Test  4  are  shown  in  Figures  7,8,  and 

9,  respectively.  In  the  shadowgraph  and  OH*  images,  a  soot  cloud  can  be  viewed  rotating  about  the  gaseous  oxygen 
core.  The  soot  clouds  corresponds  to  the  dark  regions  in  the  flow  field  that  obscure  the  view  to  the  flame.  As  mentioned 
earlier,  this  soot  cloud  surrounds  the  oxygen  core  and  blocks  the  view  to  the  flame  and  its  anchoring  regions. 
Occasionally,  highly  luminous  reacting  pockets  can  been  seen  to  attach  to  the  injector  exit  in  Figure  7.  Given  the  CH* 
and  OH*  averages  indicate  the  flame  is  separated  from  the  nozzle  but  is  visibly  attached  within  the  nozzle  in  the 
shadowgraph  average  and  instantaneous  images  in  Figure  7,  it  is  vital  that  direct  observation  of  the  flame  within  the 
injector  is  needed  to  verify  the  flame’s  anchoring  location. 

The  soot  cloud  near  the  injector  droplet  field  obscuring  the  near  field  region  of  the  flame  is  caused  by  the  design 
of  the  injector.  GCSC  injectors  are  designed  to  have  a  fuel -rich  outer  layer  that  helps  keep  the  injector  face  from 
overheating.  In  a  multi-element  configuration,  this  fuel  flow  across  the  injector  is  pushed  downstream  when  fuel  from 
neighboring  elements  meet  in  the  interelement  regimes.  This  creates  a  fuel-rich  mixture  that  breaks  down  as  it  is 
heated,  creating  a  soot  cloud  that  can  block  light  emission  in  the  UV  and  visible  spectrums. 

Observation  of  the  high-speed  videos  shows  the  soot  cloud  originates  from  the  periphery  of  the  swirling  liquid 
sheet,  which  is  consistent  with  a  fuel-rich  mixture  breakdown  being  responsible  for  the  soot  formation.  The  soot  cloud 
also  rotates  in  the  direction  of  the  swirl  and  hence  carries  the  angular  momentum  imparted  by  the  swirl.  The  generated 
soot  clouds  convect  at  much  slower  velocity  compared  to  the  oxygen  core  and  the  flame.  Unreacted  fuel  droplets  are 
also  observed  to  be  shed  from  the  periphery  of  the  liquid  sheet.  Although  droplets  were  observed  in  the  flow  field,  the 
spatial  resolution  was  not  high  enough  to  provide  droplet  statistics  for  these  studies.  The  OH*  chemiluminescence 
was  also  very  diffusive  compared  to  the  other  two  viewing  methods  and  it  is  not  fully  understood  why  this  is  the  case 
at  the  moment.  From  the  instantaneous  CH*  images,  it  was  observed  that  the  luminous  zones  were  created  in  isolated 
pockets,  indicating  that  the  flame  this  not  spatially  continuous  but  instead  reacts  in  isolated  regions.  These  isolated 
pockets  were  observed  to  extend  further  downstream  from  the  field  of  view  and  hence  no  flame  length  statistics  were 
extracted  for  these  studies. 

C.  Proper  Orthogonal  Decomposition 

Proper  orthogonal  decomposition  (POD)  was  select  to  study  the  large  structures  associated  with  the  flame  and 
rotating  soot  cloud  present  in  the  instantaneous  images.  Only  the  POD  analysis  for  Test  4  will  be  discussed  in  greater 
detail  as  the  results  for  the  other  test  conditions  yielded  similar  findings.  The  POD  analysis  decomposes  a  time  series 
of  images  into  spatially  orthogonal  structures,  with  each  spatial  mode  having  an  energy  and  temporal  component.  A 
representation  of  the  most  dynamic  flow  features  can  be  reconstructed  by  combining  the  most  dominant  POD  modes14. 
This  POD  analysis  was  performed  on  the  average  image  subtracted  fluctuating  component  (i.e.,  instantaneous  image 
minus  the  average  image)  for  each  optical  method. 

The  energy  associated  with  each  mode  from  the  shadowgraph,  CH*,  and  OH*  POD  results  are  shown  in  Figure 

10.  From  this  analysis,  there  is  no  dominant  mode  present  in  the  flow  field,  but  instead  the  modal  energies  have  a 
linear  relationship.  This  indicates  that  reconstructing  a  representative  flow  field  of  the  dynamic  processes  is  only 
possible  by  combining  hundreds  of  modes,  as  there  are  no  truly  dominant  modes.  This  lack  of  dominant  modes 
indicates  that  the  rotating  soot  cloud  is  very  chaotic  and  turbulent.  Inspecting  the  certain  spatial  modes  along  with  the 
fast  Fourier  transform  (FFT)  of  the  corresponding  temporal  mode,  the  rotation  of  the  soot  cloud  about  the  oxygen  core 
can  be  observed  along  with  its  physical  spatial  scales.  This  is  given  by  the  alternating  red  and  blue  pattern,  which 
represents  changes  in  light  intensity  caused  by  the  soot  cloud  blocking  the  view  to  the  flame.  This  pattern  is 
independently  detected  with  each  optical  viewing  method,  hence  the  soot  cloud  is  a  dominant  flow  field  structure  in 
these  studies. 
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Multiple  frequency  peaks  are  present  within  the  single  spatial  mode  shown  in  Figure  1 1 .  These  multiple  frequency 
peaks  indicate  that  the  creation  of  the  soot  clouds  and  their  rotation  about  the  oxygen  core  is  irregular.  A  similar 
observation  can  made  by  performing  an  FFT  (shown  in  Figure  12)  on  the  pixel  intensity  that  is  15  inner  injector 
diameters  away  from  the  nozzle  exit.  This  location  was  chosen  as  it  encompasses  the  average  path  taken  by  the  soot 
clouds  throughout  the  one  second  of  high-speed  images  captured.  From  Figure  12,  multiple  frequencies  are  also 
present  without  a  defined  peak,  similar  to  the  POD  temporal  modes. 

D.  Acoustically  Forced  Results 

For  the  acoustically  forced  studies,  Test  4  was  studied  in  greater  detail.  The  injector  was  forced  at  six  different 
chamber  resonant  frequencies,  corresponding  to  PN  excitation,  spanning  from  1920  Hz  to  3430  Hz.  For  each  of  the 
six  forcing  frequencies,  the  standing  wave  pressure  amplitude  was  set  to  20.68  KPa  (3  Psi).  Using  the  definition  of 
DMD  decomposition,  where 

I(x,  y,  t)  =  Real(£/j=1  Aj  exp(i27r/;At)  DMj ),  Eq.  1 

I{x,y,t)  is  the  time  series  video,  Aj  is  the  amplitude  or  modal  energy  of  component  j,  i  is  the  imaginary  number 
coefficient,^  is  the  complex  frequency  of  mode  j ,  At  is  the  time  separation  between  each  frame,  and  DMj  is  the  complex 
spatial  mode  j.  In  DMD,  each  spatial  mode  j  has  a  unique  complex  frequency  /.  We  can  define  the  sensitivity  of  the 
reacting  flow  to  the  acoustic  forcing  by  comparing  the  Aj  relative  to  its  neighboring  spatial  modes,  as  comparing  Aj 
to  its  neighbors  is  similar  to  comparing  frequency  peaks  from  power  density  plots  to  its  neighbors. 

Figure  1 3  shows  the  Aj  values  for  the  unique  frequencies  of  each  DMj ,  similar  to  a  power  density  plot.  The  results 
for  the  shadowgraph  high-speed  video  will  be  discussed  further  as  the  CH*  DMD  output  yielded  similar  results.  It 
should  be  noted  that  dynamic  mode  decomposition  was  not  applied  to  the  unforced  conditions  as  temporal  resolution 
for  these  experiments  was  not  high  enough  as  per  the  Nyquist  criteria.  For  the  forcing  cases,  however,  the  DMD 
analysis  shows  the  strongest  peaks  relative  to  the  neighboring  spatial  modes  is  strongest  for  1920  Hz  and  3060  Hz 
forcing.  The  shape  of  the  spatial  mode  can  be  seen  in  Figure  14  for  the  forcing  frequency  1920  Hz  and  3060  Hz.  For 
both  frequencies,  the  acoustic  perturbations  impose  the  same  shape  to  the  flow,  where  the  injector  experiences  in¬ 
plane  flapping  motion.  The  in-plane  flapping  is  imposed  onto  the  central  gaseous  oxygen  flow,  behaving  almost  as  a 
single  jet.  The  in-plane  flapping  motion  is  given  by  the  yellow  and  blue  bands  directly  across  from  each  other  right 
below  the  nozzle  exit.  The  effect  on  the  acoustic  perturbations  on  the  liquid  film  or  droplets  was  not  observed. 

The  in-plane  flapping  motion  can  be  observed  by  viewing  the  temporal  evolution  of  mode  j  at  the  forcing  frequency 
fj  detected  by  the  DMD  analysis.  The  temporal  evolution  of  the  spatial  mode  throughout  time  can  be  created  by  using 
Eq.  1 .  By  selecting  the  mode  j  of  interest  and  marching  At  in  time,  the  time  evolution  of  a  specific  mode  can  be  created. 
The  temporal  evolution  of  mode  j  associated  with  1920  Hz  can  be  seen  in  Figure  15.  The  spatial  mode  was  cropped 
to  focus  on  the  behavior  near  the  nozzle  exit.  The  red  and  blue  pattern  alternate  along  a  given  side  of  the  gaseous 
oxygen  core.  This  alternating  pattern  is  due  to  the  in-plane  flapping  motion  of  the  gaseous  oxygen  core.  The  physical 
interpretation  of  this  red  and  blue  pattern  is  changes  in  light  intensity  associated  with  the  flame  along  the  oxygen  core 
swaying  side  to  side. 


Summary  and  Conclusions 

Although  the  results  presented  here  cover  a  small  amount  of  the  vast  physical  processes  associated  with  GCSC 
injectors,  there  are  several  findings  from  this  experimental  study.  Among  the  findings  are  the  following: 

•  Flame  appears  to  be  separated  from  the  nozzle  exit  as  shown  in  the  CH*  and  OH* 
chemiluminescence  average  image,  but  attached  in  the  shadowgraph  images 

•  Depending  on  the  operating  J  condition,  the  flame  is  asymmetric 

•  Cone  angle  formed  by  the  flow  is  also  asymmetric 

•  Soot  cloud  is  created  at  the  outer  periphery  of  the  swirling,  liquid  sheet 

•  Soot  cloud  is  non-periodic  in  nature  and  rotates  about  the  oxygen  core,  blocking  the  flame’s 
chemiluminescence 

•  Soot  cloud  convects  slowly  compared  to  the  oxygen  core  and  the  flame 

•  From  the  POD  analysis,  the  lack  of  high  energy  modes  and  presence  of  the  multiple  strong 
frequency  peaks  indicate  that  the  soot  cloud  is  highly  stochastic  in  its  creation  and  rotation  about 
the  oxygen  core 


5 

American  Institute  of  Aeronautics  and  Astronautics 


•  The  gaseous  oxygen  core  behaves  similar  to  a  single  jet  under  acoustic  forcing 


•  From  the  DMD  analysis,  the  GCSC  injector  has  shown  a  sensitivity  to  the  frequency  of  the 
external  acoustic  field 

It  is  suspected  that  the  sensitivity  of  the  injector  to  the  acoustic  forcing  can  be  due  to  injector  geometry,  but  will  have 
to  be  investigated  further.  Future  work  will  be  identifying  if  the  momentum  flux  ratio  dominants  the  dynamics  of  the 
reacting  flow  as  it  does  for  the  cold  flow  studies.  Given  the  limited  range  of  J  values  tested  here,  no  conclusions  can 
presently  be  drawn  regarding  the  influence  the  momentum  flux  ratio  has  on  the  reacting  flow  field.  Sensitivity  and  its 
subsequent  behavior  to  acoustic  excitation  under  a  pressure  anti-node  configuration  will  also  be  investigated. 
Combined,  these  future  studies  will  extend  the  current  knowledge  of  GCSC  injectors. 
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Figure  1.  A  schematic  of  a  gas-centered  swirl-coaxial  injector, 
where  rg  represents  the  initial  gas  radius,  Lc  the  injector  cup  length, 
rp  the  gas  post  radius  at  the  end  of  the  sheltering  lip,  r0  the  outlet 
radius,  S  the  step  height  and  r  the  initial  film  thickness. 


Figure  2.  A  CAD  rending  of  the  high  pressure,  windowed  combustor. 


Figure  3.  A  schematic  of  the  optical  setup.  A  collimated  light  passed 
through  two  dichroic  mirrors  used  for  CH*  and  OH*  chemiluminescence. 
The  blue  represents  the  431nm  CH*  chemiluminescence,  the  purple 
308nm  OH*  chemiluminescence,  and  the  yellow  the  broadband 
collimated  light  for  used  for  shadowgraph. 
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Figure  4  -  (a)  False  color  CH*  average,  (b)  False  color  OH*average  and  (c)  Shadowgraph  average.  The  vertical 
green  line  is  located  at  the  nozzle  exit  plane  and  is  10  mm  in  length.  Note  that  images  are  not  at  the  same  scale. 
From  left  to  right,  the  averages  are  from  Test  1, 2,  3,  and  4,  respectively. 
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Figure  5  -  Wavelet  fused  image  of  shadowgraph  and  CH*  averages.  The  red  represents  the  CH*,  while  the  blue 
background  represents  the  shadowgraph.  The  purple  vertical  line  is  3  inner  core  diameters  from  the  nozzle 
exit.  Top  row  from  left  to  right  is  Test  1,  2  and  bottom  row,  left  to  right  is  Test  3  and4,  respectively. 
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Figure  7  -  Time  series  of  an  instantaneous  shadowgraph  images.  A  dark  soot  cloud  can  be  seen  obstructing  the 
view  to  the  center  core. 
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Figure  8  -  Time  series  of  an  instantaneous  CH*  images.  The  flame  is  longer  than  the  field  of  view.  The 
differential  pressure  transducer  is  blocking  part  of  the  view. 


Figure  9  -  Time  series  of  an  instantaneous  OH*  images.  The  flame  is  longer  than  the  field  of  view  and  the  soot 
cloud  obstructs  the  same  regions  as  the  shadowgraph. 
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Figure  10  -  POD  Energy  values  associated  with  each  POD  spatial  mode,  a)  Shadowgraph,  b)  CH*,  and  c)  OH*. 
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Figure  11  -  A  POD  Spatial  mode  and  FFT  of  the  corresponding  temporal  mode.  A)  Shadowgraph,  B) 
CH*,  and  C)  OH*.  The  alternating  pattern  is  due  to  the  soot  cloud  blocking  the  flame’s  chemiluminescence. 
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Figure  12  -  A  FFT  of  the  pixel  intensity  from  the  shadowgraph  images  taken  15  oxygen  diameters  away  from 
the  nozzle  exit. 
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Figure  13  -  The  DMD  amplitude  A  vs  the  frequency,  a)  is  1930  Hz,  b)  2380  Hz,  c)  2780  Hz,  d)  2852  Hz,  e)  3060 
Hz,  and  f)  3430  Hz. 
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a)  b) 

Figure  14  -  The  DMD  spatial  mode  extracted  from  the  high-speed  shadowgraphy  for  a)  1920  Hz,  and  b)  3060 
Hz.  The  alternating  band  is  indicative  of  a  in-plane  flapping. 


Figure  15  -  The  DMD  spatial  mode  temporal  evolution  extracted  from  the  high-speed  shadowgraphy  of  the 
1920  Hz  acoustic  perturbations.  An  in-plane  flapping  motion  is  indicative  by  the  changing  red  and  blue  pattern. 
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AF  relevant  considerations 


*  Achieving  modern  thermodynamic  efficiencies  requires 
achieving  increasingly  higher  chamber  pressures,  sometimes 
exceeding  the  critical  pressure  of  the  reactants 

—  e.g.,  liquid  rockets,  future  gas  turbines 

*  When  the  combustion  systems  are  for  propulsion,  limited 
tankage  dictates  that  on-board  propellants  be  stored  in 
condensed  form 

—  e.g.,  kerosene,  liquid  oxygen  in  rockets 

*  Combustion  systems  can  no  longer  be  designed  to  meet 
modern  requirements  without  considering  system  dynamics 

*  Combustion  dynamics  always  includes  acoustic  waves, 
which  in  enclosed  systems  can  sometimes  reach  detrimental 
amplitudes 

—  e.g.,  combustion  instabilities 
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Objectives 


•  Study  a  single  geometry  sub-scale  gas  centered,  swirl  coaxial 
injector  (GCSC)  in  a  high  pressure  environment  using  ambient 
RP-2  and  gaseous  Oxygen 

•  Study  the  reacting  flow  field  using  three  optical  methods 
simultaneously. 

*  Shadowgraph 

*  CH*  Chemiluminescence 

*  OH*  Chemiluminescence 

•  Study  the  reacting  flow  field  for  various  momentum  flux  ratios  J. 

•  Perform  a  sensitivity  study  of  the  injector  geometry  to  acoustic 
perturbations. 
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Approach 


Open  the 
loop  for 
better 
control 


Acoustic 

waves 


Thermoacoustic 
Instability 
Mechanism 

Hydrodynamic 

Unsteady  response 

Heat  release 
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Gas  Centered  Swirl  Coaxial  Injector 


Features 


*  Oxygen  rich  gas  is  injected 
through  the  center  of  the 
injector. 

*  Liquid  fuel  is  injector 
tangentially  from  the  center 
flow. 

*  A  swirling  film  is  created 
and  shrouds  the  inner 
gaseous  flow 


•  The  film  is  atomized  due  to 
shear  and  tangential  forces 
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Experimental  Facility 


Features 

-  Frequency  and  amplitude  independent  of  combustion  -  accurate  control  of  frequency  and  amp. 

-  Pressurization  independent  of  combustion  -  accurate  control  of  pressure. 

•  Subcritical  and  supercritical  pressures 

-  Precise  cryocooler  -  accurate  control  of  temperature  to  within  ±1  K. 

-  Chamber-within-a-chamber 

•  Outer  chamber  contains  pressure  -  pressure  containing  elements  remain  cool 

•  Inner  chamber  contains  acoustics  and  combustion  only  -  allows  finer  adjustment  of  inner  elements 

-  High  amplitude  piezosirens  specially  designed  for  high  pressure  J  1 

-  On-axis  windows  for  shadowgraph,  Schlieren,  chemiluminescence,  OH*  emission  Rayleigh  Index 

-  Off-axis  windows  for  PIV/PLIF  fields 

-  Fully  developed  turbulent  injector  flows  -  well  known  boundary  conditions 

-  High-speed  pressure  transducers  + - * 
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Summary  of  Forcing  Conditions 


*  Pressure  node  (PN)  and  pressure  antinode  (PAN)  at 
the  injector  location 


Imposes  transverse 
velocity  oscillation 


Imposes 

unsteady 

backpressure 


•  Forcing  frequency  ~  3000  Hz 

*  Pressure  fluctuation  amplitudes  (peak-to-peak) 
range  up  to  approximately  9  psi  (6  psi  reacting) 
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Operating  Conditions 


Study  four  J  operating  conditions  at  3.2 
MPa  (460  Psi) 

J  ~  4.4  to  6.5 

Fuel  Velocity  -  6.5  m/s 

Oxygen  Velocity  -  58.7  m/s 

Kept  fuel  massflow  the  constant  and 
varied  oxygen  to  achieve  the  desired  J 
value 


Use  the  optical  set  up  of  shadowgrpahy,  CH*  and  OH*  chemiluminescence. 


Acoustically  force  the  reacting  flow  to  at  6  resonant  pressure  node 
conditions  from  1920  Hz  to  3430  Hz.  The  pressure  node  was  kept  constant. 

Collected  1000  frames  at  1000  fps  for  average  behavior  and  then  5000  frames 
at  25kHz  for  the  acoustic  studies. 
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Unforced  Results  (J  =  4.8) 


Shadowgraph  CH*  OH* 


*  Soot  clouds  are  created  at  the  periphery  of  the  swirling  film. 
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Unforced  Results  Shadowgraph 

Average 


J  =  4.8 


J  =  5.5 


J  =  6.5 


•  The  flame  is  attached  to  the  nozzle  exit. 
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Unforced  Results  CH*  Average 


J  =  4.4  J  =  4.8  J  =  5.5  J  =  6.5 


*  A  flame  appears  to  be  separated  from  the  nozzle. 
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Unforced  Results 


J  =  4.4 


J  =4.8 


J  =  5.5 


J  =  6.5 


9  From  the  CH*  average  (in  red)  it  is  separated  from  the  nozzle  3 
inner  jet  diameters  away  from  the  nozzle. 
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Unforced  Results  OH*  Average 


J  =  4.4  J  =  4.8  J  =  5.5  J  =  6.5 


*  From  the  OH*  there  is  an  asymmetry  depending  on  J. 


DISTRIBUTION  A.  Approved  for  Public  Release;  Distribution  Unlimited.  PA  Clearance  16597  13 


Asymmetric  Cone  Angle 


J 

Left  0  ° 

Right  6  ° 

4.4 

5.27 

10.45 

4.8 

5.79 

9.94 

5.5 

6.84 

9.90 

6.5 

5.58 

7.88 

*  The  asymmetry  is  also  present  in  the  cone  angle. 
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Proper  Orthogonal  Decomposition 


In  POD  the  flow  field  is  decomposed  into  uncorrelated  coherent 
structures  with  it  corresponds  energy  content  and  temporal  evolution. 

•  Perform  proper  orthogonal 
decomposition  (POD)  on 
fluctuating  component  of  the  J=  4.8 
condition. 

•  POD  was  used  to  study  the  large 
scale  structures  associated  with 
the  reacting  flow.  Particularly  the 
soot  cloud  and  its  rotation  about 
the  central  oxygen  core. 

•  POD  was  applied  for  the  three 
optical  methods. 

•  Performed  POD  on  the  data 
collected  at  1000  Hz  over  one 
second. 


I(x,  y,t)  =  UZV 


Peter  Schmid, 201 3 
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Singular  Values 


POD  Eigen-Modes 


Shadowgraph 


*  No  dominant  modes  are  present 


CH* 


OH* 
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OH* 


POD  Eigen-Modes 


CH* 


*  No  dominant  modes  are  present 


Shadowgraph 
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Singular  Values 


POD  Eigen-Modes 


OH* 


*  No  dominant  modes  are  present 


Shadowgraph 


CH* 
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POD  Spatial  Modes 


Shadowgraph 


*  Pattern  represents  the  soot  cloud  rotation  about  the  oxygen 
core. 
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Forced  Results  (J  =  4.8) 


Shadowgraph  (25  kHz) 

Unforced  1920  HZ  3060  HZ 


DISTRIBUTION  A.  Approved  for  Public  Release;  Distribution  Unlimited.  PA  Clearance  16597  20 


Forced  Results  (J=  4.8) 


CH*  (25  kHz) 

Unforced  1920  HZ  3060  HZ 
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Dynamic  Mode  Decomposition 


Extract  spectrally-pure  temporal  modes  with  detailed  spatial  mode 
shapes 

•  Schmid  (2010)  and  Rowley  et  al.  (2009) 

•  Employ  time-averaged  amplitude  measurement  described  by 

AleniuS  (2014)  Amplitude  of  mode  at  t  =  0 

•  2000  samples  used 


n 


\ 


I(x,  y,  t)  =  Re  ^  Af  exp {Xit)Di  (x,  y) 


1=1 


Accounts  for  growth  of 
mode  in  time  as  well  as 
temporal  frequency 


Complex  spatial 
mode  shape 


Properties  of  DMD 

•  Isolates  response  of  flow  at  forcing  frequency  and  harmonics 

•  Single  modes  can  reconstruct  convective  processes  (POD  requires  two  modes) 

•  Less  efficient  at  reconstructing  signal  energy  compared  to  POD 
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Amplitude  Amplitude 


Dynamic  Mode  Decomposition 
 Shadowgraph 


1920  Hz 


2380  Hz 


2780Hz 


3430  Hz 


*  The  flow  was  receptive  to  1920  Hz  and  3060  Hz  forcing. 
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DMD  Results 


1920  Hz 


3060  Hz 


» * 


•  The  dominant  DMD  spatial  physically  represents  in-plane 
flapping  of  the  oxygen  core. 
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Conclusion 


*  Flame  appears  to  be  attached  at  the  nozzle  exit  from  the 
shadowgraph  images. 

*  Soot  clouds  block  optical  view  to  the  flame  anchoring  regions  in 
the  CH*  and  OH*  images 

*  Soot  clouds  are  created  at  the  periphery  of  the  liquid  film  and 
rotate  around  the  gaseous  oxygen  core. 

*  The  POD  analysis,  the  soot  clouds  are  large  flow  structure 
associated  with  the  flow  and  are  highly  stochastic. 

*  The  injector  has  selective  forcing  frequencies  to  which  it  is  most 
receptive. 

*  The  DMD  analysis  detects  the  in-plane  flapping  motion  of  the 
oxygen  core. 
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Window  Coating 


No  coating  Window  Coated 


DISTRIBUTION  A.  Approved  for  Public  Release;  Distribution  Unlimited.  PA  Clearance  16597  27 


